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Abstract 

In the charged pion decay, a neutrino is produced in pair with a charged lepton and they have the 
same production rate. In this paper we show that neutrinos have their own space-time correlations 
in a wide area and are detected in a different manner from charged leptons, owing to extremely 
small mass. The neutrino flux reveals a unique interference effect in the form of diffraction of 
non-stationary waves. The diffraction component of the flux shows a slow position-dependence 
and leads to an electron neutrino at short base-line regions. The electron neutrino flux at short 
distances is attributed to the neutrino diffraction and the one at long distances is to the normal 
flavor oscillation. The former depends upon the average mass-squared and the latter depends 
upon the mass-squared difference drn^. The LSND and the two neutrino experiment (TWN) 
measure and the other experiments measure (5m^. Hence they are consistent with each other. 
The neutrino diffraction would supply valuable information on the absolute neutrino mass. 
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NEUTRINO DIFFRACTION AND LSND. 



Time evolution of a charged pion state is determined by the Schrodinger equation with 
the weak interaction Hamiltonian, and the state vector of a single pion at t = becomes a 
superposition of one pion state and a state composed of a charged lepton and a neutrino at 
t = T. A decay of the pion is studied normally with a transition amplitude at the infinite T 
of plane waves, i.e., asymptotic boundary conditions. This paper presents that T-dependent 
probabilities of the decay products give valuable informations, particularly on the neutrino, 
that are not found with the asymptotic values. Since the wave function describes a time 
dependent development of the neutrino, the wave function is used to compute the time- 
dependent probability, which reveals wave phenomena. These physical quantities have been 
unable to observe experimentally with high precision before due to low statistics, but they 
are becoming possible with high statistics data of the current experiments. We study the 
neutrino detection probability in a situation where the initial pion is prepared in a plane 
wave using an S-matrix of a finite time interval. 

If the neutrino at a detector region is highly correlated, this wave is different from an 
isolated free wave and is not expressed by the plane wave of satisfying the asymptotic 
boundary condition of the S-matrix. In this region, expressing neutrinos with plane waves 



fl3, 



are not appropriate. 



and studying their detection processes with the ordinary S-matrix 
An S-matrix of a finite time interval defined by the M0ller operator at the finite time interval 
n±(T) as, S[T] = i7^T)f2+(T), is appropriate to study the finite-size correction of the 
neutrino detection probability. Since wave packets decrease rapidly at large \x\ and satisfy 
the asymptotic boundary conditions, the S-matrix of a finite time interval is expressed with 
the wave packets. One feature of the scattering of the finite time interval is that ^[T] does 
not commute with the free Hamiltonian Hq but satisfies 

[s[nHo] = ^{-^n4T)yn4T)-^n4T)^-^n4n (1) 

thus the energy is not conserved by S[T]. Final states of having the different energy from 
that of the initial state contribute to the finite-size correction of the transition probability. 
Furthermore, the finite-size correction has a universal property and is computed rigorously. 

We study the neutrino in the pion decay with 5'[T], and find that the neutrino flux has 
a large finite-size correction of a form of a diffraction. We summarize our results first. The 



fluxes of tlie neutrino and cliarged lepton in pion decays are defined with their detection 
rates and are given at a macroscopic distance in the form 



P — Pnormal + PdiffO^), (2) 

where Pnormai agrees with a normal term calculated with the standard method and Pdiff(T) 
is a diffraction term which is derived from the energy-nonconserving final states. The diffrac- 
tion components for the neutrino, P^'^jj^{T), and for the charged lepton, Pjjjj(T), are ex- 
pressed by the masses m^.,mi and the mixing matrix Ut^e, where i is the mass eigenstate, in 
the form 

p!^}j(T) = CT5^^(T,a;.J|f/,,p, P^^^ = CTg{T,uj,), (3) 

i 

2 

where C is a constant obtained later, u = and g{T,u) = g(T,ijj) — (7(00,0;), where 

Tg{T,u) = -i[ rftidt2TFTe'"^*- (4) 



^(T, u) is positive definite and decreases with a product Tu. Tui becomes large and ^(T, ui) 
vanishes at a macroscopic time for the charged leptons, but g(T,ijj^.) becomes finite for the 
neutrinos. Hence Pdiff{T) is finite in neutrinos and depends upon an average mass-squared 
''^l- Pnormal IS the normal term and includes flavor oscillations with the period determined 
by the mass-squared differences 6ml. For charged leptons the diffraction terms vanish and 
physical quantities are computed with the normal term. 

Due to the small mass, relativistic invariance, and other features, the diffraction effect 
becomes observable in the neutrinos. The detection rates become different from their pro- 
duction rates at the macroscopic distances. Especially because the energy-momentum is 
not conserved in the diffraction component, the helicity suppression mechanism does not 

work and the electron neutrino is not suppressed compared to the muon neutrino in the 
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near-detector region. This resolves the LSND anomaly 3l| and some previous experiments. 



POSITION-DEPENDENT DETECTION PROBABILITY: ONE SPECIE. 



Now we derive the diffraction term 



To show this effect being distinct from a flavor 



oscillation, the formula for one specie of neutrino is studied first. 



We suppose that a neutrino is observed through its incoherent interaction with one of 
nucleus. Then its detection amphtude in the pion decay process is expressed as, T = 
J d'^x {l,i>\Hyj{x)\7i). Here a pion is prepared at a time T^r, and a neutrino is detected 
at Tu, where the distance c(Tj, — T^r) is macroscopic. A detected neutrino is expressed 
by a wave packet that represents the nucleon wave function in a nucleus that neutrinos 

5- 



interact. The neutrino wave packet 



7| is described by the central values of the momentum 



and coordinate, (Ty,Xy) and the width, |8Nl4|. They are expressed in the form |7r) 



|pV,T^), |/,z/) = \l,'pi;h','Pi^,Xu,T^). The amplitude T is written with the hadronic V — A 
current and Dirac spinors in the form 

T = j d''xdKN^{0\Ji}^^{x)\n)u{pi)^^il - ^^HK) 

^^ipi-x+iku-{x-X^)-^{k^-p^)'^ 
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where Ni = ig {(J^/tx) ^ {mim^/ EiE^) ^ , and the time t is integrated in the region T^ < t. 
is the size of the neutrino wave packet and was estimated using the size of a nucleus. The 
Gaussian form of the wave packet is used for the sake of simplicity to obtain the finite-size 
correction in this paper. Its long-distance behavior is the same in general wave packets as 
was verified in jj]. The muon momentum of ^ p-^ — Pu and broad tail of p^ ^ p-w — Pu 
contribute to this amplitude. The former component is the normal one and the latter one 
is a diffraction component which is shown to be computable rigorously using a light-cone 
singularity of relativistic invariant systems. 

The neutrino momentum is integrated easily in Eq. and the coordinate representa- 
tion of the neutrino wave is obtained, which shows the time evolution of the neutrino wave 
function in the backward direction. At t = T,y, the wave function agrees with the Gaussian 
function of the center Xj, and at t < Tj, the position of the center is at v^{t — T^) + X^, 
which overlaps with the pion and muon wave functions. 

Integrating the space coordinates, a Gaussian function of the momenta, which shows that 

the momenta are approximately conserved, are obtained. The time is integrated in the finite 

interval T = T^, — Ttt and the amplitude is proportional to the T-dependent term, 

sin [{E{p^) + E{p^) - E{p^) - ■ % + - P^))T/2] 
E{p^) + E{p^) - E{p^) - ■ {p^ +Pf,- Ptt) 

The energy is modified by the tv dependent term and particularly the neutrino energy and 
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momentum are combined to the small value, E{p^) — Vy ■ p^ = Hence the muon energy 



can be larger than that of the energy-momentum conserved system and leads the large 
finite-size correction. 

To compute the amplitude and probability of this component rigorously, we introduce a 
correlation function and write the probability, after the spin summations are made, in the 
form 

dpi 



(27r) 



Em 



51,^12 



= ^ f d^xid''x2e-^^^^^'-^^^"Ki{5x)e'^^'''\ (7) 

Ey J 

3 — * 

where N2 = (An/a^)^ V^'^, K is a normalization volume for the initial pion, = Xj, + 
— Tj,), 6x = xi — X2, (p{6x) = py-5x and 



xe-'(P--K)'5^. (8) 

In this expression the momentum is integrated first, which is possible because the probability 
is finite and integration variables can be interchanged. 

The waves of charged lepton at the ultra-violet energy regions cause a light-cone singu- 
larity. Changing the variable io q = p-^ — pi and integrating the four dimensional g, we are 
able to extract the light-cone singularity 8{8x'^)e{6t) {4, 15| easily. The integral from the 
region g*' < leads 5{5x'^)e{5t) and the terms written by Bessel functions, while that from 
< g*^ < leads the rapidly oscillating term . The light-cone singularity is real without 
oscillation and is extended to infinity and the others either oscillate or decrease rapidly. So 
this expression is useful to find the finite-size correction of the probability which is unable 
to obtain with standard calculations of plane waves. In the first region the energy is not 
conserved, and the integral vanishes at T = cxd in fact. The second one, on the other hand, 
is that of conserving the energy and determines the quantities at T = 00. This expres- 
sion of writing the probability with the light-cone singularity converges and is valid in the 
kinematical region 2p^-p,^ < fhf, where = — mf. 

Substituting the expression of At,^i{6x) into Eq. ([7]), we have the phase factor of the 
neutrino wave at the light-cone position x = cn^{t — T,,) +X;y in the form, = E,^(t — T^) — 

— * 2 

Py\x — y^y) = ^^{t — T„) of a slow angular velocity. Next an integration on the coordinates, 
Xi and X2, and ti and ^2 in the finite T = T^, — T^, leads the slowly decreasing term g{T,u!) 



and the normal term Go- g{T,uj) is generated from the hght-cone singularity and related 
term and at T = oo, g{T,uj,y) vanishes. The normal term, Go, is from the rest. Due to the 
rapid oscillation in St, Go gets contribution from the microscopic 6t region and is constant in 
T. This term does not depend on a^, and agrees with the normal probability obtained with 
the standard method of using plane waves. In the region 2pT^-py > mf, Aj^i{6x) does not 
have a light-cone singularity and the diffraction term exists only in the kinematical region 
'^Pn'Pu < TTT'f- This region depends upon the charged- lepton mass, hence the diffraction 
terms of all three mass eigenstates converge in the union of the kinematical regions of the 
three masses, 2p^-pj, < m^. The diffraction terms is applied in this region in this paper. 
We compute the total probability next. From the integration of neutrino's coordinates 
the total volume V is obtained and cancelled with the normalization of the initial pion 
state. The total probability, then, becomes sum of the normal term Gq and the diffraction 
term ^(T, Uu), 

where N3 = STg'^di, and L = cT is the length of the decay region. P is the neutrino flux 
when it interacts with the physical state of finite size target of (7,y. At finite T, the flux has 
the diffraction component, which is caused by the superposition of waves and stable under 
the variation of the pion's momentum. At T — )■ 00, the diffraction term vanishes and the 
probability P agrees with the value of the standard calculation using plane waves. 

Now we study each term in Eq. (|9]). In the normal term, energy and momentum are 
approximately conserved, and Gq has a sharp peak at p^-Pu = ^f/'^- Hence the factor 
— 2j9^ • py in Eq. ^ becomes mf and the rate is proportional to mf. Integration of the 
neutrino's angle leads this integral independent of the angle width, as far as it include the 
narrow peak. The value is independent also from a,^, which is consistent with the condition 
for the stationary state |12|, and agrees with the value of the ordinary method, which has the 
suppression of the electron mode. On the other hand, the diffraction component is present 
in the wide kinematical region, \py\{ET, — \Ptt\) < P-n-Pv < and depends on a^,. The size 

of the nucleus of the mass number A, Oy = A^/m'^ makes the value, = 5.2/m^ for the 
^^C nucleus. Since the energy non-conserving states contribute, the diffraction component 
has a larger value of — 2pT^-py in Eq.([9]) than mf. Hence the branching ratio to the 
electron mode can become much larger than that of the normal terms and is determined by 
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FIG. 1. The neutrino energy and • pi, dependences of the fraction of Pdiff are given in 1-a 
and 1-b. The horizontal axis shows in [MeV]. ruy = 0.25 - 2.0 [eV/c^], = 2 [GeV], and 
L = 50 [m]. p-kPv dependence of the diffraction term is broad and that of the normal term is 
narrow. 

the integral of the diffraction component over the angle covered by the detector in the above 
region. This value is sensitive to the geometry of the experiment and could be much larger 
than 10"'*, if the distance between the decay region and the detector is small. 

The energy and Pn-Pu dependences of the fraction of the diffraction component are pre- 
sented in Fig.H] for = 0.25 - 2 [eV/c% E^ = 2 [GeV], and L = 50 [m]. The energy 
distribution of the diffraction term has a peak at Ey = 100 — 200 [MeV], and the Pj^-Pu 
distribution is broad in the diffraction component and is sharp in the normal component. 

From Fig. ([T]), there is an excess due to the diffraction in the short distance region and 
the maximal excess is about 0.2 of the normal term. The diffraction term is slowly varying 
with both the distance and energy. The typical length Lq of this universal behavior is 
Lo [m] = 2Eyhc/ml = 20 x Ey[50MeV]/ml[eV^ /c^]. The diffraction is stable with the 
energy, hence a measurement with a finite uncertainty of the neutrino's energy AE^, which 
is of the order 0.1 x i^j,, does not change the value much. For instance, with the uncertainty 
5 [MeV] of the energy 50 [MeV] the diffraction component does not change at all. For a 
larger value of energy uncertainty, too, the diffraction is stable under a change of the energy 
and is computed from Eq. ([9]). 
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FIG. 2. The fraction of the diffraction probabihty at a finite distance L are given for LSND 
(2-a)and TWN (2-b). The horizontal axis shows the distance in [m] and the vertical axis shows 
fractions. The two values are very different but are consistent and agree with the diffraction 
predictions, fh^ = 0.25 - 2.0 [eV/c^], = 350 [MeV/c] in LSND and 300 [MeV/c] in TWN, and 
= 60 [MeV] in LSND and 100 [MeV] in TWN. 

ELECTRON NEUTRINO IN PION DECAY. 



Formulae for general three families are easily obtained from the one species formula. 
Since the light-cone singularity in A^ i(5a;) and the diffraction term have the universal form 
that is independent of the mass of the charged lepton, and depends upon the absolute 
neutrino mass, the diffraction term to the electron neutrino is obtained from the mixing 
matrix ?7j g, Eq. ([3]). Hence the flux of the electron neutrino is written with the average mass 
ml, X]i ^(T, i^z/J|f^i,eP = g(T^,^u),^u = ^If^E^. The average value coincides with m^. if 
the 6ml are much smaller than the average value. 

Historically using neutrinos in picn decays. TWN Q proved that two neutrinos are 

different, and LSND [3| claimed observation of — )■ z/g transition with a large 6ml = 



0.1 — 1.0 [eV^/c"^]. Unfortunately this value is inconsistent with MiniBooNE jl7| which was 
designed to to test LSND and the values 6ml = 7.6 x 10~^,2.4 x 10~^ [eV^/c''] from other 
experiments jisj]. TWN found also signals of electron neutrino flux at an order of 10 — 20 
per cent, which is too large for the flavor oscillation assumption. So electron neutrino events 
of LSND and TWN are inconsistent with other neutrino experiments, if they are attributed 
to flavor oscillations. They have been puzzles for quite some time. 
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TABLE I. Geometries of experiments 



Experiment Lds Ld-reg 

TWN 30 [m] 20 [m] 

LSND 30 [m] 0.7 [m] 

CDHSW(N) 120 [m] 50 [m] 

MiniBooNE 500 [m] 25 [m] (diameter = 12.2 [m]) 

CDHSW(F) 600 [m] 50 [m] 

KARMEN - Stopped pion or muon decay 

TABLE II. Experimental values and theoretical values 

Experiment PvjPv^, (Exp) Pdiff / Pnormai (Th) 

TWN 0.18 0.17 

LSND (2.6 ± 1.0 ± 0.5) X 10-3 2.8x10^3 

CDHSW(N) ? 0.2-0.5 

CDHSW(F) ? 0.02 - 0.05 

MiniBooNE, KARMEN ^ 10"^ 



Now we compare the theoretical value of electron neutrino events with experiments. 
The normal component is negligibly small due to helicity suppression and the diffraction 
component Pdiff is finite and is studied. The neutrino flux in the decay pipe area is computed 
from the diffraction formula and the neutrino expressed by the wave packet of determined 
by the detector propagates freely in the next area between the decay pipe and the detector 
region. The distribution is spread with p^-Pv and has different geometrical behaviors from 
the normal component in the detector. Effects of geometries are important. The geometries 
of TWN and LSND are similar and the lengths between the pion source and the neutrino 
detector, Ld-s, and those of the decay region, Ld-reg, are given in Table HI 

0.7 [m] of Ld-reg is the length of the decay region in air, while the whole decay region 
is 1.8 [m]. Since the length between the detectors and the sources of MiniBooNE and 
CDHSW(F) Q is as large as 10 times those of TWN and LSND, the contributions of 
diffraction terms are ^ of TWN and LSND. At CDHSW(N), the length is 4 times of TWN 
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Distance between target and detector [mj 

FIG. 3. Electron neutrino flux of the neutrino diff'raction is compared with flavor osciflations. The 
horizontal axis shows the distance in [m] and the vertical axis shows the fraction of the electron 
neutrino. The diffraction agrees with the LSND. The value from Am,^ of the standard oscillation 
parameters becomes negligible at 30 M. To fit with the experiment by the flavor oscillation, a 
sterile neutrino of Am^ ~ XieVjc?^'^ is necessary. 



and LSND but the energy is higher and about the same magnitude of the diffraction term 
is obtained. The diffraction term comes from the tree diagram and is suppressed strongly 
in matter, KARMEN [20I]. The diffraction component is important in decay in air but is 
complicated in decay in matter. Its magnitude will be studied in a future publication. We 
compute the fraction of the diffraction components by taking into account the geometry. In 
LSND, the decay region in air is about 0.7 [m] long and the detector is located 30 [m] away 
from the decay region. In this case, the ratio becomes 0.28 percent at = 1 [eV/c^]. In 
TWN, the decay region is about 20 [m] long and the detector is located 10 [m] away from 
the decay region. The neutrino detection angle is much wider and the ratio becomes about 
17 percent at = 1 [eV/c^]. In CDHSW, the energies of pion and neutrino are higher and 
iron is used for the detector. Including these effects, theoretical values are obtained. 

The ratios Pdiff /Pnormai of our theory and Pu^/Pu^ of experiments are summarized in 
Table [771 Uncertainties in TWN are not known and ignored. The fractions of the diffraction 
components obtained in this way are presented in Fig. [2] for the mass of the neutrino, rhi, = 
0.25 - 2 [eV/c^], and the pion momentum = 350 [MeV/c] in LSND and 300 [MeV/c] in 
TWN, and the neutrino energy = 60 [MeV] in LSND and = 100 [MeV] in TWN. The 
fractions of the theory agree with those of the experiments using the same rrij^. 
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The magnitude of electron neutrino flux agreed with LSND in flight and TWN but not to 
others including MiniBooNE. Since matter effect in LSND at rest is complicated, we study 
only LSND in flight here. In Fig.3, the electron neutrino fraction from the neutrino diffrac- 
tion is compared with those of the flavor oscillations of standard parameters of Am^ and 
of a sterile neutrino. The neutrino diffraction explains naturally the excess of the electron 
neutrino. Thus the puzzles are resolved and consistent understanding of the neutrino exper- 
iments are obtained. Furthermore, information on the absolute neutrino mass is obtained 




SUMMARY AND IMPLICATIONS. 

It is found that the neutrino detection rate has the unusual flnite-size correction in the 
macroscopic distance. Instead of the standard S-matrix ^[oo], which is useless to study the 
flnite-size correction, we studied the S-matrix of the flnite time interval S[T] and computed 
the flnite-size correction of the detection probabilities which are measured in experiments. 
The asymptotic boundary conditions must be satisfled for the scattering or decay processes, 
hence the wave packets are necessary for ^[T] and are used. One feature is that S[T] does 
not commute with the free Hamiltonian Hq and both flnal states of conserving the energy 
and those of non-conserving energy contribute. The rate from the former is independent 
of both T and and agrees with that of the standard method, while that from the latter 
depends on T and Uj, and does not satisfy various relations of the former such as the helicity 
suppression, symmetry between the neutrinos and charged leptons, and others. 

The electron neutrino flux due to the diffraction was computed in the geometries of 
LSND and TWN and excellent agreements between the theory and experiments are obtained. 
Because the electron neutrino flux due to diffraction is determined by the average neutrino 
mass m^, the excesses of the electron neutrino suggest that rh^ is around 0.25 — 1.0 [eV/c^]. 
If the rhi, is in this range and of the experiments are used, LSND and TWN electron 
events are attributed to the diffraction component but not to the flavor oscillation. Thus the 
controversy related with LSND and TWN are resolved. Some information on the absolute 
neutrino mass is obtained, furthermore. 

The conflrmation of the neutrino diffraction would be made in a new experiment from the 
unique behaviors of the neutrino spectrum on the Pt^'Pu, energy, and angle. Since the helicity 
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suppression does not work in the diffraction component, the flux of the electron neutrino 
becomes about 10—20 % of that of the muon neutrino at the maximum value. Furthermore it 
has the unique properties that are very different from the normal term. If clear distinctions 
are observed on the dependences of these variables and agreements with the theory are 
found, the neutrino diffraction will be confirmed. Since at much larger distance than the 
above length, the diffraction component disappears and only the normal flavor oscillation 
terms remain, the experiment to measure the excess due to the diffraction component should 
be made at a short distance. 

The present results may exhibit violations of two important relations, the unitarity and 
lepton number conservation. However that is not the case and they are not violated. First, 
from the behavior of the diffraction component that the detection rate is decreasing with 
T, the neutrino number looks decreasing and an unitarity is violated superflcially. Since the 
neutrino is detected at the distant position, the detection time is delayed. The retarded and 
interference effects made the detection probability deviate from the neutrino number. Hence 
the decreasing behavior of -P(T) is in-connected with both the decreasing neutrino number 
and unitarity. In the text, P(T) is computed with S[T] that satisfies the proper boundary 
condition of the experiments and the unitarity condition 5'[T]^5'[T] = 1. If the T dependent 
neutrino probability were computed with <S'[oo], this would not have decreased with T. 
However, the S[oo] does not satisfy the proper boundary condition for the neutrino at T 
and is neither applied nor compared with the probability measured by the experiment. S[T], 
on the other hand, satisfies the proper boundary condition and is used for comparisons with 
the experiments at the finite T. One of the important findings of the present work is that 
S[T] gives the different results from S[oo] even in the macroscopic T. Second, asymmetric 
behaviors of g{T,uj) lead that the neutrino detection rates are very different from those of 
the charged leptons. This is not in-consistent with the lepton number conservations but 
means that the detection probability is different from the neutrino number owing to the 
finite-size correction caused by retarded and interference effects. When either the neutrinos 
or the charged leptons are observed, their detection rates are different. Nevertheless, when 
both of neutrino and charged lepton are observed simultaneously, they have the same rates, 
due to the lepton number conservations. 

The shorter life time of the pion in the situation when the neutrino is observed at a 
small distance than that of the free pion has not been tested and will be confirmed in a 
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future experiment. Then if the muon is measured simultaneously, this muon shows the same 
behavior from the lepton number conservation. The detection rate of the muon depends on 
the boundary condition on the neutrino. 

In this paper we ignored the higher order effects such as the pion life time and the pion 
mean free path in studying the quantum effects. We will study these problems, the muon 
decay, and other large scale physical phenomena of low energy neutrinos in subsequent 
papers. 
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